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We construct quantum circuits for measuring the commuting set of vertex and plaquette operators 
that appear in the Levin- Wen model for doubled Fibonacci anyons. Such measurements can be 
viewed as syndrome measurements for the quantum error-correcting code defined by the ground 
states of this model (the Fibonacci code). We quantify the complexity of these circuits with gate 
counts using different universal gate sets and find these measurements become significantly easier 
to perform if n-qubit Toffoli gates with n — 3, 4 and 5 can be carried out directly In addition 
to measurement circuits, we construct simplified quantum circuits requiring only a few qubits that 
can be used to verify that certain self-consistency conditions, including the pentagon equation, are 
satisfied by the Fibonacci code. 



I. INTRODUCTION 

The ground states of certain two-dimensional lattice 
Hamiltonians of a type first introduced by Kitaev^can be 
used as quantum error-correcting codes known as surface 
codes. Quantum information can be stored and protected 
using these codes when they are defined on lattices with 
holes (defects) P Fault-tolerant gates can then be carried 
out either transversally or by deforming the code in order 
to braid these defects while staying entirely within the 
code subspacepHS One downside to using the Kitaev sur- 
face codes, for which defects behave as Abelian anyons, 
is that to realize a universal set of fault-tolerant gates 
at least one gate using a resource costly "magic state" 
distillation process 6 is required. The same is true for 
fault-tolerant quantum computation using the so-called 
color codesP"^ Nevertheless, quantum computation us- 
ing these surface codes has a number of appealing fea- 
tures, notably the need for only nearest-neighbor gates 
between qubits in a two-dimensional array and high erro r 
thresholds, e.g. ~ 1% for the Kitaev surface code.CSML^ 

Recently Konig, Kuperberg, and Reichardt 13 (KKR) 
outlined a method for fault-tolerant quantum computa- 
tion using non- Abelian surface codes. These codes, which 
are defined mathematically in terms of the Turaev-Viro 
topological invariants for 3-manifolds, 14 can be viewed 
physically as ground states of Levin- Wen models,^ two- 
dimensional lattice models which generalize the Kitaev 
model.^ These models can be used to realize so-called 
"doubled" versions of any consistent anyon theory, in- 
cluding theories of non- Abelian anyons for which braid- 
ing is universal for quantum computation. The simplest 
such universal anyons are the Fibonacci anyons. Here 
we refer to the corresponding Levin- Wen model as the 
Fibonacci Levin- Wen model and, following KKRp^ refer 
to the ground states of this model as the Fibonacci code. 
As shown in Ref. 13, when using the Fibonacci code, 
Fibonacci anyons can be associated with holes in the lat- 
tice subject to certain boundary conditions and proper 
initialization. These Fibonacci anyons can then be used 



to encode logical qubits and universal quantum compu- 
tation can be carried out purely by braiding themp"^ 
without the need for magic state distillation. 

The Levin- Wen models are defined by a set of com- 
muting vertex and plaquette projection operators which 
act on qubits (more generally, qudits) associated with the 
edges of a two-dimensional trivalent lattice. When using 
the ground states of these models as quantum codes it 
will be necessary to continually measure these vertex and 
plaquette operators in order to check for errors, which 
would then have to be corrected without disturbing the 
quantum information stored in the topological degrees 
of freedom of the code. For the Kitaev surface code, 
quantum circuits which can be used to measure these op- 
erators are known and are fairly straightforward.^ For 
either an n-sided plaquette, or a vertex where n edges 
meet, these measurement circuits each require a single 
initialized syndrome qubit which is measured after car- 
rying out n controlled-NOT (CNOT) gates. 

The simplicity of the quantum circuits used to measure 
the vertex and plaquette operators for the Kitaev surface 
code reflects the Abelian nature of this code. It is nat- 
ural to ask how complex the quantum circuits need to 
be to measure the vertex and plaquette operators for the 
non-Abelian Fibonacci code. In this paper we present 
explicit quantum circuits for performing such measure- 
ments. These circuits are built in part out of smaller cir- 
cuits which carry out unitary transformations which have 
been described both in KKRp^ and, in the context of en- 
tanglement renormalization, in Ref. |2l] Our goal here is 
to explicitly construct these circuits in terms of standard 
elements (Toffoli gates, CNOT gates and single-qubit ro- 
tations) in an attempt to quantify their complexity. 

The purpose of this work is not to argue that non- 
Abelian surface codes are viable competitors to the Ki- 
taev surface code. Indeed, we share the view of many 
in the field that quantum computation using the Kitaev 
surface code, given its clear advantages over other fault- 
tolerant quantum computation schemes, may well pro- 
vide the best prac tical route to building a functioning 
quantum computer Here our goal is the more modest 
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FIG. 1: (Color online) (a) Example of a trivalent lattice (in this case a honeycomb lattice) on which the Levin- Wen model 
can be defined. For the Fibonacci Levin- Wen model a qubit is associated with each edge. A particular state which satisfies 
the vertex constraint Q v = 1 on each vertex for this model is shown. Thick edges indicate qubits in the state |1), thin edges 
indicate qubits in the state |0). (b) Action of the vertex operator Q v on the three qubits on the edges connected to a trivalent 
vertex, and of the plaquette operators on the 2n qubits associated with an n-sided plaquette. 



one of making a first pass at determining the complexity 
of syndrome extraction for the significantly less well un- 
derstood Fibonacci code, which we believe is of intrinsic 
interest in its own right. An additional goal of this work 
is to begin developing a "dictionary" for translating the 
mathematical structures which appear in general anyon 
theories into interesting quantum circuits, some of which 
we find require only a few qubits and might feasibly be 
carried out experimentally in the near future. 



II. LEVIN- WEN MODELS AND THE 
FIBONACCI CODE 

The Levin- Wen models^ are defined on two- 
dimensional trivalent lattices such as the hexagonal lat- 
tice shown in Fig. [IJa). The degrees of freedom of the 
models are associated with lattice edges which can take 
on a finite number of labels. These labels can, in gen- 
eral, be oriented, meaning for each label i there is a dual 
label i*. If i = i* then the edge is unoriented. For the 
Fibonacci Levin- Wen model there are only two labels 
and 1 and the edges are unoriented (0 = 0*, 1 = 1*). 
Thus, for this model, as for the Kitaev surface code, we 
simply associate a qubit with each edge of the lattice. 
The two states of each qubit |0) and |1) then correspond 
to the two labels and 1, respectively. 

For a given trivalent lattice the Levin- Wen Hamilto- 
nian has the form 



# = -£Qv-$> p 



(i) 



Here Q v and B p are projection operators associated with 
the vertices (labeled v) and plaquettes (labeled p) of the 
lattice. 

The vertex operator Q v acts on the three qubits associ- 
ated with the edges connected to vertex v and is diagonal 
in the standard {|0), |1)} basis. (Here we focus on the Fi- 
bonacci Levin- Wen model and so only consider the case 
when a single qubit is assigned to each edge.) If these 



qubits are in the states \j) and \k) the result of ap- 
plying Q v is determined by the tensor 5^ (see Fig.[TJb)) 
which, for the Fibonacci Levin- Wen model, is given by, 



fi if 

\ ot 



ijk = 000,011,101,110,111 
otherwise. 



(2) 



The plaquette operator B p is significantly more com- 
plex than Q v . For example, for a hexagonal plaquette, 
B p acts on the six qubits on the edges of plaquette p 
in a way determined by the state of the six qubits on 
the edges connected to the plaquette. B p is therefore a 
twelve-qubit interaction (in general a 2n-qubit interac- 
tion for an n-sided plaquette). For the Fibonacci Levin- 
Wen model the precise form of the plaquette projection 
operator is, 



(*2 + ^J)> 



(3) 



where B p for s 



and 1 are plaquette operators asso- 
ciated with the label s and <p — (\/5 + l)/2 is the golden 
ratio. The action of B p on an n-sided plaquette is shown 
in Fig. [ljb) where, 
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Here the six-indexed tensor F^ , along with forms 
the basic data of a so-called tensor category — the math- 
ematical framework for a general anyon theory, in this 
case the theory of Fibonacci anyons. The F and S tensors 
satisfy certain self-consistency conditions which, among 
other things, guarantee that the operators B p and Q v 
all commute with each other! 15 * 23 ! Note that since the Fi- 
bonacci Levin- Wen model is unoriented, in ([5| we have 
assumed i = i* for all labels. The precise form of the F 



tensor for this model is given in Sec. IV 



When using the ground states of the Levin- Wen model 
as quantum error-correcting codes the commuting vertex 
and plaquette projection operators Q v and B p should be 
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viewed as stabilizers. The code space is then defined by 
the requirement that Q v = 1 on each vertex and B p = 1 
on each plaquette. For the Fibonacci code the constraint 
Q v = 1 projects the Hilbert space onto the space spanned 
by states in which edges in the state |1) form branching 
loop configurations (see Fig. []Ja)), while the plaquette 
constraint B p = 1 leads to particular quantum superposi- 
tions of these states. As described in KKR,^when these 
code states are defined on lattices with holes that have 
certain boundary conditions on their edges, these holes 
(or defects) can realize a "doubled" version of the anyon 
theory characterized by the F and 6 tensors. For the 
Fibonacci code this means that these defects can encode 
two types of Fibonacci anyons with opposite chiralities. 
As further shown in KKR, 13 with proper initialization 
these defects can be forced to encode Fibonacci anyons 
of a particular chirality. These chiral anyons can then be 
used to encode qubits and braided in order to carry out 
universal quantum computation. 

In this paper we focus on the problem of how to mea- 
sure the stabilizers Q v and B p for the Fibonacci code. In 
the passive approach to the Levin- Wen model envisioned 
in KKRp^ rather than engineering the Levin- Wen Hamil- 
tonian to realize the Fibonacci code it will be necessary 
to continually measure these operators in order to detect 
errors which can then be corrected. 



III. QUANTUM CIRCUIT TO MEASURE Q v 

The measurement of Q v for the Fibonacci code is 
straightforward and not significantly more difficult to 
carry out than the analogous measurement for the Ki- 
taev surface code. A quantum circuit which carries out 
a quantum non-demolition measurement of Q v is shown 
in Fig. [2] The circuit acts on the three qubits associated 
with a given vertex as well as a fourth syndrome qubit 
initialized in the state |0). After carrying out the circuit 
the syndrome qubit is measured. If it is found to be in 
the state |0) then Q v = 1 for this vertex and the vertex 
constraint is satisfied, if not then Q v = and the vertex 
constraint is violated. 

The most difficult part of the Q v circuit to carry out 
is likely to be the four-qubit Toffoli gate which performs 
a NOT gate on the syndrome qubit if and only if the 
state of each of the three vertex qubits is |1). (Here and 
throughout it should be understood that an n-qubit Tof- 
foli gate is a gate with n — 1 control qubits and one target 
qubit.) This four-qubit Toffoli gate is the first of several 
n-qubit Toffoli gates required in our constructions, all of 
which are directly related to the non-Abelian nature of 
the Fibonacci code. Here this gate is needed to allow for 
the loop branching associated with the fact that 5m = 1. 

In what follows we will be interested in quantifying 
the complexity of the quantum circuits we construct. Of 
course the notion of quantum circuit complexity is some- 
what ill-defined and depends, among other things, on 
what we take as our primitive gate set. This in turn 
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FIG. 2: (Color online) Quantum circuit which can be used to 
measure Q w for the Fibonacci code. 



will depend on the particular hardware of the quantum 
computer being considered. 

Accurate three-qubit Toffoli-class gates have recently 
been been carried out experimentally using supercon- 
ducting qubit d 24 " 26 ! and trapped ions. 27 Motivated by 
this, we take one primitive gate set to consist of three- 
qubit Toffoli gates, CNOT gates and single-qubit rota- 
tions. An n-qubit Toffoli gate can then be carried out 
using 4n — 12 three-qubit Toffoli gates if n — 3 additional 
qubits are available™ These additional qubits need not 
be initialized and their states are left unchanged once the 
full n-qubit Toffoli gate is carried out. Thus nearby code 
qubits which are not being acted on directly by the opera- 
tor under measurement can be used. With this construc- 
tion we can count the total number of three-qubit Toffoli 
gates (or, simply, Toffoli gates), CNOT gates and single- 
qubit rotations required to carry out a given circuit. For 
the case of the four-qubit Toffoli gate appearing in our 
Q v circuit this count gives 4 Toffoli gates. The total gate 
count for our Q v circuit is then 4 Toffoli gates and 3 
CNOT gates. This can be contrasted with the analogous 
circuit for the Kitaev surface code which, when acting on 
a trivalent vertex, would require only 3 CNOT gates (it 
is, in fact, identical to the circuit shown in Fig. [2] with 
the four-qubit Toffoli gate removed)™ 

For a second gate count we assume that the n-qubit 
Toffoli gates which appear in our circuits are themselves 
primitive gates. By this count, our Q v circuit consists 
of 1 four-qubit Toffoli gate and 3 CNOT gates. We 
note that there are proposals for carrying out single-step 
n-qubit Toffoli-class gates using trapped ions™ super- 
conducting qubits,^ and neutral atoms interacting with 
cavity photons™ in addition, it has been observed that 
these gates are efficiently achieved if one of the qubits 
has n available quantum levels™ Of course n-qubit Tof- 
foli gates can also be simulated using the usual primi- 
tive gate set consisting of CNOT gates and single-qubit 
rotations.^ However, as we have seen with our Q v mea- 
surement circuit, and as will become more clear in what 
follows, the ability to directly carry out accurate n-qubit 
Toffoli gates (with n = 3, 4 and 5) will give a strong ad- 
vantage when carrying out quantum computation using 
the Fibonacci code. 

Despite requiring a four-qubit Toffoli gate, the Q v 
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FIG. 3: (Color online) (a) An F-move, a five-qubit unitary 
operation denned in terms of the tensor F"$*, . (b) Action of 
an F-move on the abstract trivalent lattice of the Fibonacci 
code which illustrates the decoupling of this lattice from the 
physical qubits. In this example the qubits (open circles) are 
arranged in a Kagome lattice and lie on the edges of an initial 
trivalent (hexagonal) lattice. After the F-move the edges of 
the new trivalent lattice must be distorted if they are forced 
to coincide with the physical qubit lattice. 

measurement circuit shown in Fig. [2] is relatively simple, 
reflecting the simplicity of the vertex operator. In what 
follows we turn to the more difficult problem of mea- 
suring the plaquette operator B p . For this case a brute 
force approach to constructing a circuit which measures 
the appropriate operator acting on the edges of a pla- 
quette for each possible state of the edges connected to 
that plaquette is problematic. Fortunately, there is a use- 
ful resource which simplifies the problem greatly — the 
F-move. 



IV. F-MOVE 

When using the Fibonacci code, the physical qubits 
of a quantum computer may be fixed in space and may 
even form a rigid lattice. However, this physical lattice 
need not be the same as that formed by the edges of the 
abstract trivalent lattice used to define the code. Indeed, 
as emphasized in KKR,^ this abstract trivalent lattice 
should be thought of as fluid and constantly changing 
throughout the computation. These changes are accom- 
plished by carrying out F- moves, processes which locally 
redraw the trivalent lattice while reassigning the physical 
qubits to new lattice edges and carrying out an appro- 
priate unitary operation. 

Specifically, when carrying out an F-move five edges 
of the lattice are redrawn as shown in Fig. j3^a) while 
a unitary transformation determined by the six indexed 
tensor F^, (the same F tensor which appears in ^) 
is applied to the five qubits associated with these edges. 
This five-qubit unitary is a controlled operation on the 
qubit labeled e in Fig.[3ja) contingent on the states of the 
other four qubits (labeled abed). The usefulness of the 
F-move here derives from the fact that if one starts in a 



ground state of a given Levin- Wen model on a particular 
trivalent lattice then, after performing an F-move, the 
resulting state will be a ground state of the new Levin- 
Wen model denned on the new trivalent lattice. 21 This 
is true even though this lattice has decoupled from the 
physical qubits, as illustrated in Fig. j3^b). 

It was shown in KKR 13 that the ability to decouple 
the abstract trivalent lattice from the physical qubits 
with F-moves is an important resource for carrying out 
quantum computation using the Fibonacci code. For ex- 
ample, by carrying out sequences of F-moves one can 
deform the code to perform Dehn twists on the trivalent 
lattice which can then be used to braid defects encod- 
ing Fibonacci anyons.^ Since the braiding of Fibonacci 
anyons is universal for quantum computation, this means 
that one can perform a universal set of gates while stay- 
ing inside the Fibonacci code subspace without the need 
for magic state distillation. 

The F-move for the Fibonacci code is represented 
graphically in Fig. [4] This figure, together with Fig.^a), 
can serve as a definition of the F tensor for Fibonacci 
anyons. The effect of carrying out an F-move is only 
shown for those states which satisfy the vertex constraint 
(i.e. for which Q v = 1 for all vertices). When defining 
the Levin- Wen models, the F tensor is assumed to van- 
ish when acting on those states which violate the vertex 
constraint. 15 Here we will assume before applying any F- 
move that it has been verified that Q v = 1 on each rele- 
vant vertex of the initial trivalent lattice. The structure 
of the F-move then guarantees that the vertex constraint 
will continue to be satisfied on the new trivalent lattice. 

A quantum circuit which acts on five qubits at a time 
and which carries out the F-move defined in Fig. [4] for 
those states satisfying the vertex constraint is shown in 
Fig. |5| In this fi sure the labels abode refer to the same 
labels shown in Fig.|3|a). Although it is not immediately 
apparent from its structure, one can readily check that 
this circuit has the symmetries of the F tensor^ (e.g., 

F cdt> = Kbt' = F dct>)- Note also that the circuit squares 
to 1 (since F 2 = 1, see below), so the same circuit can be 
used for the inverse transformation. As described above, 
this F circuit carries out a particular operation on the 
qubit labeled e depending on the state of the other four 
qubits labeled abed which are themselves left unchanged 
at the end of the circuit. The F circuit can therefore be 
viewed as a generalized Toffoli-class gate. Because the 
four control qubits are not equivalent, it is important to 
label these qubits in our F circuit as we have done in the 
green box in Fig. [4] This notation will be useful when we 
embed the F circuit into larger circuits acting on more 
than five qubits. 

At the heart of the F circuit is the five-qubit 
controlled-F gate where F is the 2x2 unitary matrix 
acting on qubit e when a = b = e = d = 1, 




The remaining Toffoli gate and CNOT gates take care of 
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FIG. 4: (Color online) F-move for Fibonacci anyons. Under this F-move a unitary transformation is performed on the qubit 
associated with the edge which goes from horizontal to vertical conditioned on the state of the qubits on the other four edges. 
As in Fig. [I] thick lines indicate edges in the state |1) and thin lines indicate edges in the state |0). Only those states which 
satisfy the Q v = 1 constraint are shown. 



all other cases for which the outcome is essentially fixed 
by the vertex constraint. As stated above, this circuit 
is designed to carry out an F-move only on those states 
which satisfy the Q v = 1 constraint on all vertices. In 
what follows we will always assume it has been verified 
that the vertex constraint is satisfied before applying the 
F circuit. 

Figure [5^b) shows how to carry out the five-qubit 
controlled- F gate using a five-qubit Toffoli gate and two 
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single-qubit rotations. This simple construction is pos- 
sible because F 2 = 1 and detF = — 1. As for the four- 
qubit Toffoli gate appearing in the measurement circuit 
for Q v , the appearance of this five-qubit Toffoli gate 
can be traced back to the fact that loops are allowed 
to branch in the Fibonacci code and is a direct conse- 
quence of the non-Abelian nature of this code. Using the 
construction of Ref. [28] described above this five-qubit 
Toffoli gate can be carried out using 8 conventional Tof- 
foli gates. The total gate count for our F circuit is then 9 
Toffoli gates, 4 CNOT gates and 2 single-qubit rotations. 
Alternatively, if we treat n-qubit Toffoli gates as primi- 
tives, the gate count is 1 five-qubit Toffoli gate, 1 Toffoli 
gate, 4 CNOT gates and 2 single-qubit rotations. Given 
the importance of carrying out F-moves when using the 
Fibonacci code pa the ability to accurately carry out this 
five-qubit Toffoli gate can be viewed as an important ex- 
perimental threshold for realizing this type of quantum 
computation. 
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FIG. 5: (Color online) (a) Quantum circuit which carries out 
an F-move for the Fibonacci code (the 2x2 matrix F is 
given in Eq. pn. The labels abode refer to the same labels 
in Fig. |3ja). (b) Five-qubit controlled-F gate expressed in 
terms of a five-qubit Toffoli gate. Here R(±0y) = e ±i6(J y /2 are 
single-qubit rotations about the y axis with — tan -1 
for which R(0y)XR(-0y) = F. 



V. PENTAGON EQUATION 

The F-move satisfies an important self-consistency 
condition known as the pentagon equation. The pen- 
tagon equation can be represented as a sequence of F- 
moves on a seven-edged trivalent lattice as shown in 
Fig. |6ja). In a quantum computer, the lattice edges 
would be associated with qubits, labeled 1 through 7 in 
Fig. [6ja). As one follows this sequence of F- moves, the 
trivalent lattice is repeatedly redrawn while the qubits, 
which can be considered fixed in physical space, are reas- 
signed to the new lattice edges after each F-move. By the 
time one has gone all the way around the pentagon the 
trivalent lattice has returned to its original form. How- 
ever, the qubits associated with two of the edges (labeled 
5 and 6 in the figure) are swapped. 

The process of carrying out this sequence of five F- 
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a) 1 2 3 4 b) 




FIG. 6: (Color online) (a) The pentagon equation, a self-consistency condition which the F-move must satisfy. As shown 
here, the pentagon equation corresponds to a series of F- moves which take a particular 7 edged lattice (upper left) back to 
an identical lattice (lower left) while two of the qubits associated with the lattice edges are swapped. Here and in subsequent 
figures the edges associated with the initial state before each F-move are color coded as in Fig. [3] (b) The pentagon equation 
as a quantum circuit identity. The sequence of F-moves shown in (a) are carried out by repeatedly applying the F circuit 
defined in Fig. [5] The labels abode in each green box refer to the labels in Fig. [5] The circuit equality holds provided the 
vertex constraint Q v = 1 is satisfied on all three vertices in the initial lattice. In the figure, the triplets of numbers given below 
"Qv = 1" in the red box indicate the qubits which meet at these vertices. 



moves and the resulting qubit swap can be translated 
into the quantum circuit identity shown in Fig.[6^b). We 
refer to the left-hand side of this identity as the pentagon 
circuit. The solid green rectangles in the pentagon circuit 
represent the five-qubit F circuit shown in Fig. [5] and the 
corresponding abode labels are the same as the labels 
shown in Fig. [5] Again we assume that before carrying 
out the pentagon circuit it has been verified that Q w = 1 
on each of the two vertices of the initial trivalent lattice. 
It is only for this case that the circuit identity shown 
in Fig. [6^b) holds (for clarity these vertices are labeled 
by their associated qubits inside the red box under the 
equals sign in this figure) . 

In the pentagon circuit two of the qubits (qubits 5 and 
6) are acted on while the remaining qubits play the role 
of control qubits. Simpler quantum circuits can be con- 
structed by fixing these five effective control qubits to be 
in a particular state. For example, if we fix all the qubits 
except for 5 and 6 to be in the state |1) then the pentagon 
circuit reduces to the simple two-qubit circuit shown on 
the left-hand side of the circuit identity in Fig. (7) This 
simplified pentagon circuit consists of five controlled- F 



gates with alternating control qubits, and the net effect 
of this sequence of gates is a SWAP gate. Note that when 
qubits 5 and 6 are both in the state |0) and all other 
qubits are in the state |1) the vertex constraint is vio- 
lated in the full seven-qubit pentagon circuit. However, 
in this case the simplified pentagon circuit merely carries 
out the identity operation, which is consistent with swap- 
ping the two qubits. Therefore the expression shown in 
Fig. [7] is an exact circuit identity, regardless of the vertex 
constraint. 

We note the resemblance of this circuit identity t o the 
familiar three CNOT construction of the SWAP gatePEU 
In our case, the circuit identity shown in Fig. [7] repre- 
sents the nontrivial part of the pentagon equation which 
uniquely fixes the form of the matrix F (up to an ar- 
bitrary and irrelevant phase choice for the off-diagonal 
matrix elements). We envision that this circuit identity 
may be useful for calibrating the F operation. For ex- 
ample, one can imagine tuning F until it can be verified 
by quantum process tomography that five controlled- F 
gates with alternating control qubits indeed produce a 
SWAP gate. 
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FIG. 7: (Color online) Simple two-qubit circuit identity ob- 
tained by setting the five effective control qubits (qubits 
1,2,3,4, and 7) in the pentagon circuit identity shown in 
Fig.^b) to the state 



VI. QUANTUM CIRCUIT TO MEASURE B p 

We now turn to constructing a quantum circuit to 
measure the plaquette operator B p . To do this we use 
a method inspired by the entanglement renormalization 
scheme of Ref. [2TJ The essential idea is that through 
a sequence of F-moves any n-sided plaquette can be re- 
duced to a 1-sided plaquette with a single external line, 
i.e. a "tadpole." One such sequence of F-moves which 
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FIG. 8: (Color online) Reduction of a hexagonal plaquette to a tadpole through a sequence of F-moves. 



reduces a hexagonal plaquette to a tadpole is shown in 
Fig. [8j Note that the final F-move in this sequence acts 
on four qubits rather than five. A quantum circuit which 
carries out this reduced F-move, obtained by identifying 
the qubits labeled a and d in the circuit shown in Fig. |5j 
is shown in Fig. [9j (Gate counts for this reduced F cir- 
cuit: 5 Toffoli gates, 4 CNOT gates, and 2 single-qubit 
rotations, or 1 four-qubit Toffoli gate, 1 Toffoli gate, 4 
CNOT gates, and 2 single-qubit rotations.) 

It was shown in Ref.|2T]that the plaquette operator B p 
commutes with F-moves, i.e. after each F-move shown 
in Fig.[8]the value of B p is unchanged even as the plaque- 
tte size is reduced. This is equivalent to the statement 
that if we start with a plaquette in a ground state of 
the Levin- Wen model (meaning Q v = 1 on each vertex 
and B p = 1 for the plaquette) then, after each F-move, 
the qubits will continue to be in the ground state of the 
Levin- Wen model for the new lattice. Thus, after each F- 
move, it will still be true that Q v = 1 on each vertex and 
B p = 1 on the reduced plaquette. This means that after 
performing the "disentangling" reduction of the n-sided 
plaquette to a tadpole one need only measure B p for the 
tadpole to measure B p for the original plaquette. Since 
the tadpole only consists of two qubits this measurement 
is straightforward. 

For a tadpole there is a unique eigenstate of B p with 
eigenvalue lj^ 



Here the first qubit is the external line (tail of the tad- 
pole) and the second qubit is the 1-sided plaquette (head 
of the tadpole). The two-dimensional Hilbert space of 
states orthogonal to \ip Bp=1 ) which satisfy the vertex con- 
straint will then have B p = 0. This space is spanned by 
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|^£? P =o,a) 
IV^p=o,6) 

To measure B p for this simple two-qubit system we 
first rotate the head qubit of the tadpole so that it is in 
the state |0) if B p = 1 and in the state |1) if B p = 0. 
This can be done by carrying out a single-qubit rotation 
S on the head qubit if and only if the state of the tail 
qubit is |0) where^ 



S = 



(9) 



This transformation corresponds to the diagram shown 
in Fig. lCTa) and is defined in terms of the tensor S% b , 
which is equal to the matrix S when a = and for which 
Si! = 1 (the case Sl b , with b = or b' = violates 
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FIG. 9: (Color online) Reduced four-qubit F-move obtained 
by identifying the qubits labeled a and d in Fig. [5] 
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FIG. 10: (Color online) (a) S transformation acting on a two- 
qubit tadpole. The tensor S% h , is defined in the text, (b) 
S circuit which carries out an S transformation. The 2x2 
matrix S is given in Eq. [9] Here R(±py) = e ±ip(7y/2 are 
single-qubit rotations about the y axis with p = tan -1 _1 
for which R(py)XR(—py) = S. (c) Quantum circuit which 
uses the S circuit to measure B p for a two-qubit tadpole. 
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FIG. 11: (Color online) Quantum circuit which can be used to measure B p for the Fibonacci code on a hexagonal plaquette 
based on the plaquette reduction shown in Fig. [8] It must be verified that Q v = 1 on each of the six vertices of the plaquette 
before carrying out the circuit. 



the vertex constraint). A quantum circuit which carries 
out this transformation (and its inverse since the circuit 
squares to 1) is shown in Fig. [TO^b). This circuit can be 
carried out with 1 CNOT gate and 2 single-qubit rota- 
tions. Like the F circuit, this simple construction is pos- 
sible because S 2 = 1 and det^ = — 1. If the two tadpole 
qubits are initially in the state \i/jb p =i) the result of car- 
rying out this circuit is the state |0) |0). If the two tadpole 
qubits are initially in the two-dimensional Hilbert space 
spanned by the states {I^Bp^o? a), ^^=0? b)} then after 
carrying out this circuit they will be in the space spanned 
by the states {|0)|1), |1)|1)}. In either case the state of 
the second qubit, i.e. the rotated head of the tadpole, 
will be equal to 1 — B p . 

After carrying out the S circuit on the tadpole, a 
CNOT gate can be done with the head qubit as the con- 
trol qubit and a syndrome qubit, initialized to the state 
|0), as the target qubit. The syndrome qubit can then 



We again emphasize that the circuit shown in Fig. 11 
only measures B p correctly if the vertex constraint Q v = 
1 is satisfied on each vertex of the initial plaquette at the 
start of the circuit. If the vertex constraint is violated on 
any of these vertices then by definition B p = for the 
plaquette^ but the circuit will, in some cases, give the 
wrong result of B p = 1. To see this, consider the action 
of this circuit on the full 2 2n -dimensional Hilbert space 
of the 2n qubits associated with an n-sided plaquette, in- 
cluding those states which violate the vertex constraint. 
From the structure of the circuit, which performs a uni- 
tary transformation on 2n qubits and then measures the 
state of a single qubit to determine B p , it is clear that the 
dimensionalities of the Hilbert spaces for which B p = 1 
or ^ = would both be 2 2n_1 , i.e. half that of the 
full Hilbert space. However, once the vertex constraint is 
taken into account the Hilbert space is greatly reduced. 
The dimensionalities of the projected Hilbert spaces for 



be measured and if the result is then B„ = 1 for the whicn Qv = 1 on each of the n vertices and B p = 1 or 



tadpole (and hence for the original plaquette), and if the 
result is 1 then B p = 0. 

After measuring B p for the tadpole, the final step is 
to reconstruct the full plaquette. This can be done by 
undoing the S circuit on the tadpole and then undoing 
the F- moves. Putting everything together the resulting 
measurement circuit for the case of a hexagonal plaquette 



is the palindromic circuit shown in Fig. 11 In this circuit 
the notation is the same as in the pentagon circuit, with 
each box corresponding to either the full or reduced F 
circuit, or the S circuit, and the letters labeling the var- 
ious "inputs" as defined in Figs. |5|9[ and IlO} From the 



structure of the circuit it is clear how this construction 
generalizes to the case of an arbitrary n-sided plaquette. 



B p = for the plaquette are Dim[£? p = 1] = i^n-i an d 
D\m[B p = 0] = i^n+i, respectively, where F n is the nth 
Fibonacci number (F = 0, Fi = 1, F 2 = 1, F 3 = 2, 
etc.). For the case of a hexagonal plaquette this means 
the full 4096 = 2 12 dimensional Hilbert space of twelve 
qubits is projected down to a space of dimensionality 
322 = F11 + F13 = 89 + 233 with an 89 dimensional space 
of states satisfying the plaquette constraint with B p = 1. 
The reader will be reassured to know we have numeri- 
cally checked that the circuit shown in Fig. [TT] performs 
the correct measurement of B p on this projected space. 

It should be noted that the requirement that Q v = 1 
on each vertex before measuring B p may cause prob- 
lems when extracting error syndromes. For example, if a 
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FIG. 12: (Color online) (a) Sequence of F-moves which pulls a tadpole through a line, (b) Four-qubit quantum circuit which 
initializes a tadpole with an S circuit, carries out the sequence of two F-moves shown in (a), and then performs another S 
circuit so that measuring qubit 3 would yield B p for the new tadpole. The tadpole is initialized to a state with B p = 1 or 
depending on whether the initial state of qubit 4 is |0) or |1), respectively. The circuit equality holds provided Q v = 1 on 
the vertices of the initial lattice. As in Fig. |6jb) these vertices are labeled inside the red box. The 2x2 matrix U is given by 
Eq.IlOlin the text. 



faulty measurement of Q v gives 1 for a particular vertex 
on a plaquette, but the actual value of Q v is for that 
vertex, then, as described above, the B p measurement 
circuit for the plaquette will, in some cases, give B p = 1 
even though the correct value (as it is for any plaquette 
in which a vertex constraint is violated) is B p = 0. In 
this paper we have not addressed the important question 
of whether it is possible to extract error syndromes for 
the Fibonacci code fault tolerantly, nor the question of 
precisely how these errors would be corrected. Our goal 
has been to construct circuits which, in the absence of 
errors, can be used to measure Q v and B p in order to 
begin to get a measure of their complexity. 

We can now give our final gate counts for measuring 
B p . If we choose to reduce all n-qubit Toffoli gates to 
conventional three-qubit Toffoli gates (using 4n — 12 Tof- 
foli gates, following Ref . [28l as described in Sec. |III| ) then 
we find that our procedure for an n-sided plaquette (with 
n > 2) requires 18n-26 Toffoli gates, 8n-5 CNOT gates 
and 4n single-qubit rotations. Alternatively, if we con- 
sider n-qubit Toffoli gates as primitives, then our proce- 
dure requires 2n — 4 five-qubit Toffoli gates, 2 four-qubit 
Toffoli gates, 2n - 2 Toffoli gates, 8n - 5 CNOT gates 
and 4n single-qubit rotations.^ Not surprisingly, this is 
significantly more demanding than the analogous require- 
ment for the Kitaev surface code, for which only n CNOT 
gates are needed to measure the plaquette operator for 
an n-sided plaquette. 

Finally we note that there are, of course, many differ- 
ent ways to reduce a given plaquette to a tadpole using 
F- moves, all of which can be used to measure B p and 
some of which will be more "parallelizable" than others. 



VII. A SIMPLE EXAMPLE 

One of the motivations of the present work is to find 
simple quantum circuits which might feasibly be carried 
out in the near term and which begin to test some of the 



key properties of the Fibonacci code. We have already 
seen one example of such a circuit, the sequence of five 
controlled- F gates which results in a two-qubit SWAP 
gate discussed in Sec. [V| This circuit is a simplified ver- 
sion of the full seven-qubit pentagon circuit shown in 
Fig. [6^b) and can potentially be used to calibrate the F 
operation. In this section we give a similar example — a 
four-qubit circuit which first initializes a tadpole into a 
state with either B p = 1 or 0, and then pulls this tadpole 
through a line using F-moves to produce a new tadpole 
which can be measured to verify that the value of B p has 
not changed. As for the pentagon circuit, this four-qubit 
circuit can be simplified to a two-qubit circuit which, in 
this case, can be used to calibrate the S operation. 

The sequence of operations we consider is illustrated in 
Fig. [l2^a). The system consists of a four-edged trivalent 
lattice and so uses four qubit s, labeled 1 through 4 in 
the figure. Initially two qubits (1 and 2) are assigned 
to edges which form a line and the other two qubits (3 
and 4) form a tadpole attached to this line. As always, 
in what follows we assume that it has been verified that 
Q v = 1 on each of the two vertices of this lattice at the 
start of the process. 

The first step is to initialize the tadpole in a state with 
either B p = 1 or 0. Then, using two F-moves, as shown 
in Fig. [l2^a), the tadpole can be pulled through the line. 
The F- moves preserve B pi and so the intermediate state 
of this process is a 2-sided plaquette which has been ini- 
tialized either into the code space if B p = 1 or outside 
of the code space if B p = 0. After the tadpole has been 
pulled through the line, the two qubits forming the initial 
tadpole have swapped places — the head of the tadpole 
is now the tail and vice versa. If B p is now measured for 
the new tadpole the result should yield the same value of 
B p that the tadpole was initialized to at the start of the 
process. 

The left-hand side of the circuit identity shown in Fig- 
ure 12 'b) is a four-qubit circuit which carries out the 



procedure described above. If qubit 4 is initially in the 
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FIG. 13: (Color online) (a) Simplified two-qubit circuit iden- 
tity obtained by setting qubits 1 and 2 to the state |1) in the 

'b). (b) Equivalent circuit 
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circuit identity shown in Fig. 
identity obtained by moving the two NOT gates from the left 
side of the identity shown in (a) to the right side. 



state |1 — a) then the first S circuit initializes the tad- 
pole in a state with B p = a. A reduced F circuit then 
carries out the first F-move and produces a 2-sided pla- 
quette with B p = a. Next, a second reduced F circuit 
carries out the second F-move producing a new tadpole 
with B p = a but with the head and tail of the tadpole 
interchanged. Finally, after carrying out an S circuit on 
this tadpole the state of qubit 3 will be 1 1 — a) . 

Note that if qubit 4 is initially in the state |1) so that 
the tadpole is initialized to a state with B p = then 
qubit 3 can initially be in either the state |0) or |1) while 
still satisfying the vertex constraint. After the first S cir- 
cuit on the left-hand side of Fig. [l2]fo) is carried out the 
tadpole will then be placed in a quantum superposition 
of | -0^=0? a) and |^£ p =o,fr) (See Eqns. (FF| and (fsl) ) . At 
the end of the circuit, after being pulled through the line 
formed by qubits 1 and 2, the tadpole will still be in the 
two-dimensional B p — Hilbert space, but the particular 
superposition will in general have changed. Direct calcu- 
lation shows that if qubits 1 and 2 are both in the state 1 1) 
then the operation acting on the two-dimensional B p = 
space when pulling the tadpole through the line is given 

by 



U = 



(10) 



Otherwise, if either qubit 1 or qubit 2 (or both) are in 
the state |0) the state of the final tadpole will be the 
same as that of the initial tadpole with the head and tail 
qubits swapped. These cases are all accounted for by the 
SWAP gate and four-qubit controlled- U operation on the 
right-hand side of the circuit identity in Fig. [l2|b). As 
for the pentagon circuit, this identity only holds when 
Q v = 1 on the two vertices of the initial lattice (again 
these vertices are labeled inside the red box in the figure) . 

This four-qubit circuit, which essentially represents ini- 
tializing a 2-sided plaquette into a state with a given 
value of B p and then producing a state which can be 
measured to determine B p after carrying out a different 
F-move than the one used to initialize it, is much sim- 
pler than the full circuit for measuring B p for a hexagonal 



plaquette. However, it still involves the four-qubit Toffoli 
gate which appear in the reduced F circuit. As for the 
pentagon circuit, a simpler two-qubit circuit identity can 
be found by fixing the states of the qubits which act effec- 
tively as control qubits (qubits 1 and 2 in Fig. 12). If we 
fix these qubits to both be in the state |1) we obtain the 
simplified two-qubit circuit identity shown in Fig. 13 'a). 

This circuit identity can be simplified further by mul- 
tiplying both sides on the left and right by NOT gates 
which act on the top and bottom qubits, respectively, to 
obtain the equivalent circuit identity shown in Fig. 13 'b). 
Note that if the initial state for the circuit shown in 
Fig. [l3^a) is |1)|Q), where the first qubit is the top qubit 
(qubit 3 in Fig. 12), then the vertex constraint is not 



satisfied for the full four-qubit circuit. However, the sim- 
plified circuit identity is readily seen to be satisfied in this 
case. For all other cases the vertex constraint is satisfied, 
and so it follows that the expression shown in Fig. [l3^a) 
and the equivalent expression in Fig. 13 [b) are exact cir- 
cuit identities, independent of whether or not the vertex 
constraint is satisfied. 

The key action of the two-qubit circuit on the left-hand 
side of Fig. [H|b) occurs when the tadpole is initialized 
in a state with B p = 1 for which the tail qubit must 
start in the state |0). For this case, after pulling the 
tadpole through the line the new tadpole must again be 
in the state with B p = 1. Thus, after accounting for the 
removal of the two NOT gates, this circuit must take the 
state |1)|0) to the state |0)|1). 

Like the five controlled-F SWAP circuit in Fig. [7| 
which can be used to calibrate the F matrix, the circuit 
identity shown in Fig. [l3^b) can be used to calibrate the 
S matrix. Once F has been fixed by the pentagon circuit, 
the requirement that the circuit on the left-hand side of 
Fig, 0jb) takes the state |1)|0) to the state |0)|1) fixes 
the form of the matrix S (up to an overall phase which is 
irrelevant for our purposes). Note that in performing this 
calibration it is not necessary to carry out a full quan- 
tum process tomography. It is sufficient to verify that 
the circuit identity holds for the initial state |1)|0). For 
this case, only the SWAP gate on the right-hand side is 
relevant since the controlled-X£7 X gate enters only when 
the initial state of the second qubit is |1). 



VIII. CONCLUSIONS 

In this paper we have constructed explicit quantum cir- 
cuits for measuring the vertex and plaquette operators, 
Q v and B pi in the Fibonacci Levin- Wen model. These 
operators can be viewed as stabilizers for the Fibonacci 
codep^a surface code for which defects can behave as Fi- 
bonacci anyons — the simplest non-Abelian anyons for 
which braiding alone is universal for quantum compu- 
tation. While the Q w measurement is not significantly 
more difficult than the analogous measurement for the 
Kitaev surface codes (for which the defects behave as 
Abelian anyons), the B p measurement scheme we present 
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here is significantly more difficult than its Abelian coun- 
terpart. While the present scheme is almost certainly 
not the most efficient one for performing this measure- 
ment, given the complexity of the operator B p it is likely 
that even the most efficient schemes will require a large 
number of primitive gate operations. This cost in circuit 
complexity will then need to be weighed against the gain 
of not requiring magic state distillation. The situation 
is somewhat analogous to comparing the relative mer- 
its of performing topological quantum computation with 
Ising anyons (which requires magic state distillation) to 
Fibonacci anyons 

It is clear that further work will be needed before such 
a direct comparison of the resources needed to carry out 
fault-tolerant quantum computation using the Fibonacci 
code with that using the Kitaev surface code will be 
possible. While recent progress strongly suggests that 
the Kitaev surface code is the most promising from the 
practical point of view of trying to build an actual fault- 
tolerant quantum computer, we believe it is too early to 
rule out the possibility that the Fibonacci code may have 
some practical implications. Even if it does not, we be- 
lieve the Fibonacci code is of intrinsic interest, in part 
because computing with it can be viewed as essentially 
simulating a non- Abelian state of matter on a quantum 
computer. 

Our measurement circuit for B p is built out of circuits 
which realize F-moves and the action of the S matrix 
on a trivalent lattice. In addition to our measurement 
circuits we have also given simpler circuits built out of 
these F and S circuits. The first is a seven-qubit circuit 



which can be used to verify that the F circuit satisfies the 
pentagon equation, as well as a simpler two-qubit circuit 
which contains the nontrivial content of this equation and 
fixes the form of the F matrix. The second is a four-qubit 
circuit which uses the S circuit to initialize a tadpole (1- 
sided plaquette) in a state with either B p = 1 or 0, carries 
out a sequence of F-moves to pull the tadpole through a 
line, and then produces a state which can be measured 
to determine B p for the new tadpole. For this circuit we 
have also given a simpler two-qubit circuit which, once 
F has been fixed by the pentagon circuit, fixes the form 
of the S matrix. These simple two-qubit circuits (Figs. [7] 
and 13) may be useful for calibrating the F and S oper- 



ations. 

A recurring theme in this work has been the need for 
n-qubit Toffoli gates (with n = 3, 4 and 5) when com- 
puting with the Fibonacci code. These n-qubit Toffoli 
gates arise as a natural consequence of the non- Abelian 
nature of this code. We believe the possibility of using 
non-Abelian surface codes such as the Fibonacci code 
for fault-tolerant quantum computation provides further 
motivation for developing experimental techniques to di- 
rectly carry out accurate n-qubit Toffoli-class gates. 
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